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Abstract We analyzed the feedback of AGN jets on 
cooling flow clusters using three-dimensional AMR hy- 
drodynamic simulations. We studied the interaction 
of the jet with the intracluster medium and creation of 
low X-ray emission cavities (Bubbles) in cluster plasma. 
The distribution of energy input by the jet into the sys- 
tem was quantified in its different forms, i.e. internal, 
kinetic and potential. We find that the energy associ- 
ated with the bubbles, {pV + JP^/il ~ I))? accounts 
for less than 10 percent of the jet energy. 

Keywords Cooling flow, galaxy cluster, AGN jet, 
bubble, feedback: general 



1 Introduction 

During the last two decades clusters with cooling flows 
have been extensively studied. Prior to results obtained 
by XMM-Newton and Chandra it was widely believed 
that in these systems up to lOOOM© yr~^ of gas was 
moving from the hot phase T > 10^ K to low temper- 
atures. In this picture extensive emission by Fe XVII 
ions was predicted, in contradiction with spectra ob- 
tained by XMM-Newton (iPeterson et all 120031). The 



reasons for this conflict between theory and observa- 
tion are discussed in papers by Omma et .al (2004), 
Binney (2005) and Pizzalato & Soker (2005). 

It is now widely accepted that feedback from AGN 
jets is responsible for slowing or halting the cooling 
of gas in galaxy groups and clusters (Churazov et al. 



2000[|Birzan et al.||2004l [Binney et al.,2007| ), and it is 
likely that this phenomenon also takes place in ellip- 
tical galaxies at smaller (spatial and temporal) scales 
( Binney fc TabQrl|T995l ) . 
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AGN jets reheat hot gas by doing mechanical work 
during the inflation of low-density bubbles. The Chan- 
dra observatory has now detected bubbles as regions of 



low X-ray emission in a large number of systems ( Carilli 



et al. 1994^ David et al. 2001^ McNamara et al. 2000 



Fabian et al. 



2001). Most, but not all bubbles are as- 



sociated with non-thermal radio emission ( Birzan et al 



2004), and when both radio- loud and radio-quiet bub- 



bles occur, the radio-quiet bubbles are further from the 
center. Thus it appears that bubbles are radio-loud 
for a period after inflation, and a combination of adi- 
abatic expansion of the plasma and radiative losses by 
the most energetic electrons makes the bubble radio- 
quiet. 

In some cases a synchrotron-emitting jet runs from 
the galactic nucleus to the bubbles, so jets are definitely 
involved in bubble formation. It is unclear to what 
extent bubbles are driven by the relativistic cores of the 
jets rather than the slower and more massive flows that 
probably surround the cores. We are also uncertain 
whether the jets are driven by the black hole's accretion 



disk or emerge from its ergosphere (Blandford & Znajek 



1977), but the observation of fast jets in microquasars 
that are known to be powered by neutron stars suggests 



that an ergosphere does not need to be involved (Nipoti 
erar]|2QQ5l ). 



2 Method 

To investigate how an AGN jet interacts with cluster 
gas, we performed a series of three-dimensional hydro- 



dynamic simulations using the code Enzo (Bryan 



Norman 1997). Because of the inhomogeneity of the 
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plasma, we used the adaptive grid refinement (AMR) 
technique to subdivide the simulation domain, which 
was a cube L = 614 kpc on a side. Our coarsest grid 



2 



has 16 cells in each direction, and in the l^^ level of re- 
finement the effective number of grid points is increased 
to 16 X 2^, where / = 0, . . . , 6. The cluster core is at 
the 6*^ refinement level, so the grid would have 1024^ 
cells if the entire domain were refined to this level. We 



used the piecewise parabolic method PPM (Woodward 



& Collela 1984) Riemann solver to advance the hydro- 



dynamic equations. 

In previous simulations the effects of jets have been 
taken into account by one of three techniques: (i) an 
inner boundary is introduced on which infiow bound- 



ary conditions can be set (Basson & Alexander 2003); 
(ii) spherical cavities are carved out of the ambient 
medium and filled with much less dense material that 
provides a similar pressure dChurazov eFaL]|2QQl] ) ; (iii) 
at some arbitrarily chosen locations, energy is added 
at some specified rate (Briiggen fc Kaiser 2002 ). Our 
technique differs from approach (iii) in that we add 
mass, momentum and energy on two small, centrally 
located "injection disks" of radius rjet = 3kpc (Omma 
et al-pOOl ). The ratio of the rates of injection of mass, 
rhjet = 1^0 yr~^, and momentum specify the jet speed 
Vjet = 4 X lO^kms" , which is sub-relativistic but su- 
personic. As the simulation develops, the region in 
which the injected energy is largely thermalized moves 
away from the center, just as the hot spot of a real sys- 
tem should. Hence we believe our technique is more 
realistic than either technique (ii) or technique (iii), 
and it avoids the introduction of the unphysical inner 
boundary required by technique (i). We implement this 
scheme by solving the equations 



dp 
di 



(1) 



- V • (pl?(8)^) = -Vp-pV$ + x^jet^jet i (2) 



described by the NFW model (Navarro et al. 1995). 



The mass within a sphere of radius r is 

r/vs 



M{r) = ^nScPcrt 



ln(l + r/rs) 



1 + r/r 



(5) 



where pc is the critical density of the universe at the 
cluster redshift, is the scale radius of the cluster and 
6c is the overdensity parameter ( David et al.||2001 ) 



Sc = 



3IP 
8^ 
200 



ln(l + r2ooAs) 



^200 

l+r20o/^^s 



(6) 
(7) 



where r2oo is the radius within which the average den- 
sity is equal to 200pc- In our simulations Sc and were 
equal to 7.48 x 10^ and 77kpc, respectively. The initial 
electron density was 



ne(r) = 0.1 1 



-2.7 



cm 



(8) 



50kpc 

This provides a good fit to the Hydra cluster density 



profile reported by David et al. (2001). As initial con 



dition we assume that the system is in hydrostatic equi- 
librium and deduce the initial temperature profile from 
the assumption that the plasma is an ideal gas. 

To have a self-consistent model for jet power, one 
should assume that an AGN's outbursts are conditioned 
by the accretion rate of the black hole lurking at the 
center of the cluster. The accretion rate will depend on 
the rate of infiow of matter toward the cluster core due 
to radiative cooling loses, which form a self- sustained 
loop dBinney fc Tabor|1995[ |Cattaneo fc Teyssier|2007] ). 
Here we assumed fixed power and duration of the jet 
based on observational estimates of known AGN jets 
(Table 1). 
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-V-p^-p^- V^ + X^jet-Jrad (3) 



where p, p, and E are the mass density, velocity, 
thermal pressure and total energy, respectively, ^ is the 
gravitational potential, £^jet = \rnv'^^^^ jVad is rate per 
unit volume of radiative energy loss, and x is a window 
function that vanishes outside the injection disk. On 
the disk 



-3r^/2r 



r < r 



jet 







r > rjet, 



(4) 



where r is distance from the jet axis. The cluster bary- 
onic matter is assumed to evolve in a dark-matter halo. 



3 Results 



The cluster is allowed to cool passively for 3 x 10^ yr 



using line and continuum radiation ( [Sutherland fc Do- 
pita |1993[ ) . The cooling function A is represented in 
(Fig. [l| as a function of the gas temperature for half 
solar metalicity. 



Table 1 Parameters of AGN jet 



Jet 


life span [Myr] 


Power [erg/s] 


1 


25 


1Q45 


2 


2 X 25 


5 X 10^^ 
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Fig. 1. — Cooling function as function of gas tempera- 
ture at ^Zq. 



Fig. 4. — The volume of the bubble as function of the 
density threshold. 



The rate of radiation of energy by a shell of radius r 
and width dr is: 



jrad(r,t) = Fsheii(r,r + dr)n^(r,t)A(Te(r,t)), 



(9) 



where Fgheib and A are the shell volume, electron 
number density and cooling function (Fig. [T]), respec- 
tively. Under the assumption that the plasma is opti- 
cally thin, this energy is removed from the computa- 
tional domain. In the absence of AGN feedback, the 
system suffers a cooling catastrophe after 3.7 x 10^ yr. 

Fig. |2] shows He and Tg averaged over spherical shells 
of radius r as functions of r and t. The left panel is for a 
simulation in which there was only one outburst, while 
the right column is for a simulation with two outbursts 
lasting 25 Myr each, and separated by a quiescent pe- 
riod of 25 Myr (Table 1). 

The upper row of Fig.[2]shows that the density profile 
increases to its maximum just before the jet ignites, be- 
cause radiative losses cause the central plasma to cool, 
contract and draw plasma in from outlying regions. The 
jet immediately creates a low-density bubble (blue color 
near the cluster core). This bubble then moves forward, 
driven by both jet momentum and buoyancy (Fig[3|. 
The bubble volume is determined by locating cells in 
which the density is smaller than the averaged density 
at that radius by a factor /. Fig. [4] shows the resulting 
volume as a function of /, and on the basis of this plot 
we set / = 0.3. 

The temperature profiles in the lower row (Fig. 2) 
show strong heating during the jet phase, and this heat- 
ing can be seen to offset the cooling process for about 50 
Myr. Results not presented here show that short, pow- 
erful jets heat the gas within 10 kpc of the cluster center 



less effectively than a longer-lived lower-power jet, re- 
gardless of whether the latter has a single long spell of 
activity (~ 50 Myr) or has two 25-Myr blasts separated 
by a 25-Myr rest period. The reason is that powerful 
jets (FR II type) tend to deposit their energy far from 
the cluster core dOmma et ai;][2004l ). Of the cases pre- 
sented in Fig. 2, the jet with two cycles couples more 
efficiently with the gas within the inner cluster region. 
Although much of the gas uplifted during the first cycle 
falls back during the 25 Myr jet quiescence, the bulk of 
the energy injected in the second cycle passes to large 
radii along the channels that the jet carved through the 
dense inner ICM during its first cycle (Fig. 3). 



4 Conclusion 

From the measured bubble volumes and the pressure 
in the bubbles, we found the energy associated with 
the bubbles, /7bubbie = pV + 7P^/(7 ~ 1)^ where pV 
is the work done by the jet in creating the bubble and 
7pV/(7 — 1) is the thermal energy within the bubble, 
and 7 = 5/3 appropriate to a non-relativistic plasma. 
We found that /7bubbie is about 6 times less than the 
energy input from the jet (Binney et al. 2007| . Most 
of the energy not associated with bubbles is divided 
between the kinetic and potential energy of the gas that 
is either displaced by the bubbles or entrained in their 
wake. Some additional energy has been thermalized in 
the shocks that move ahead of the bubbles. The kinetic 
and potential energy that surrounds the bubbles will 
eventually be converted into heat through dissipation 
of the long-lived gravity waves in which this energy is 
invested. 
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One cycle jet 



Two cycles jet 
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Fig. 2. — Map of radial temperature and density profiles as function of time for one cycle jet (left column) and 
two cycles jet (right column). 




Fig. 3. — Density map along the xy plane {x is the jet axis) and bubble evolution, one cycle jet (left column); 
two cycles jet (right column). The evolution of the bubble (dark blue region) at different times 25, 50 and 73Myr 
after jet ignition. 
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